Ferrodistortive instability at the (001) surface of half-metallic manganites 
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We present the structure of the fully relaxed (001) surface of the half-metallic manganite 
Lao.7Sro.3Mn03, calculated using density functional theory within the generalized gradient ap- 
proximation (GGA). Two relevant ferroelastic order parameters are identified and characterized: 
The tilting of the oxygen octahedra, which is present in the bulk phase, oscillates and decreases 
towards the surface, and an additional ferrodistortive Mn off-centering, triggered by the surface, de- 
cays monotonically into the bulk. The narrow d-like energy band that is characteristic of unrelaxed 
manganite surfaces is shifted down in energy by these structural distortions, retaining its uppermost 
layer localization. The magnitude of the zero-temperature magnetization is unchanged from its bulk 
value, but the effective spin-spin interactions are reduced at the surface. 

PACS numbers: 68., 73.20.-r, 75.47.Lx, 75.70.-i 
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A good understanding of the surface physical proper- 
ties in colossal magneto-resistant hole-doped manganese 
oxides is highly desirable for future applications of these 
promising materials in magnetoresistive devices or spin- 
tronics. [l| The metallic phase of L&i- X A^MnOa (A =Sr, 
Ca), obtained for x ~ 0.3 is of particular interest be- 
cause of its high spin polarization at the Fermi level. 
There is some controversy, however, on whether this op- 
timally doped phase is a half-metal, with a gap for one 
spin direction, and thus with fully spin-polarized carriers. 
Early spin-resolved photoemission experiments observed 
such ideal polarization Q, although electron tunneling 
experiments show traces of partial occupation of minor- 
ity spin states Since these techniques probe the bulk 
structure through the surface or interface, (as do most 
techniques for bulk studies with variable degrees of sur- 
face sensitivity), the origin of these discrepancies could 
lie not only in the difficulty in preparing smooth surfaces 
at an atomic level, but also in the fact that the proper- 
ties of the surface may markedly deviate from those of 
the bulk. The primary goal of this work is to determine 
the structural changes induced by surface termination at 
the (100) surface of Lai_ x Sr x Mn0 3 (LSMO). In addi- 
tion, we discuss the influence of these structural changes 
on the electronic properties. 

The theoretical description of the half-metallic man- 
ganites is challenging given the strong correlations of the 
transition metal d electrons, as well as the complex chem- 
istry and large range of possible structural distortions 
from the high symmetry perovskite structure. For ex- 
ample, ab initio studies based on the local spin density 
approximation (LSDA) give poor agreement with the ex- 
perimental lattice constant; as a result many first prin- 
ciples calculations to date have been performed for cu- 



bic structures with the experimental pseudocubic lattice 
constant and without structural relaxations. However, 
the strong electron-lattice interactions are known to play 
an important role in the electronic properties, and a full 
description of the lattice distortions, particularly at sur- 
faces or interfaces where local distortions might be more 
pronounced, is desirable. Furthermore, LSDA fails to ob- 
tain a half-metallic state for x — 0.3 4j, and it has been 
argued that a more accurate treatment of electron cor- 
relation, such as LSDA+C/ or self-interaction corrected 
(SIC) is required, particularly for the charge-ordered or 
insulating phases @, [f| . 

Here we present ab initio results for a fully relaxed 
LSMO (001) surface obtained using density-functional 
theory within the generalized gradient approximation 
(GGA). Our choice of exchange-correlation functional 
was motivated by our recent GGA calculations for bulk 
Lao.7Sr .3Mn03 (LSMO), 0] in which we found good 
agreement between the theoretically optimized pseudocu- 
bic and experimental lattice constants, and qualitatively 
correct half-metallic electronic behavior. Since our focus 
is on structural, rather than on electronic properties, and 
since we are interested in the metallic phase of LSMO, 
which reduces the errors related to the non-locality of 
the exchange-correlation hole, we elect not to include 
higher order descriptions of correlation effects such as in 
the GGA+U method. Indeed, it was recently shown Q 
that the use of GGA+U does not significantly change the 
structural properties. However the limitations of today's 
GGAs for treating highly-correlated systems should be 
kept in mind, in particular when we discuss changes in 
the surface electronic spectral properties based on Kohn- 
Sham eigenvalues [§]. 

Our calculations were performed using the Siesta im- 
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plementation of density functional theory 13, II | ■ Core 
electrons were replaced by norm-conserving pseudopo- 
tentials [l2j], and finite-support numerical atomic or- 
bitals were used to describe the valence electrons, at 
the double-C polarized level Semicore states for 

Ti (3s 2 and 3p 6 ), Sr (4s 2 and %> 6 ), and La (5s 2 and 
5p 6 ) were included in the valence. The spin-polarized 
gradient-based PBE exchange-correlation functional was 
used [3] , with a grid cutoff of 320 Ry for integrations in 
real space. Atomic positions were relaxed until the forces 
were smaller than 25 me v/A. 

To validate our method in obtaining relaxed atomic 
structures, we first studied the properties of bulk LSMO. 
We built a cell that is a y/2 x \[2 x 2 repetition of the 
ideal perovskite cubic cell, and performed a variable cell 
relaxation of the crystal structure. Our calculated lattice 
parameters are within 2 % of the experimentally reported 
values, and our Mn-0 bond lengths (1.96 A) are even 
closer to experiment (1.95 A) [l5| . The tilt angle (6.6°) 
that we obtain for the optimized structure is in excel- 
lent agreement with the experimental value (6.8° [HI). 
Since these distortions are believed to be the source of 
the strong electron-phonon coupling in manganites [Til ] 
(they lift the degeneracy of e g levels, as recently shown 
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in photoemission experiments 17j |. and reduce the elec- 
tronic energy for single occupancy of this state) their 
accurate representation in our calculations is essential. 

To study the surface, we consider a supercell structure 
with 4x2 repetitions of the original cubic cell in the sur- 
face afe-plane, eleven atomic layers (5.5 unit cells) along 
the direction perpendicular to the surface, and Mn02 ter- 
mination at both sides of the slab [Hj]. We begin with a 
pseudocubic structure and allow for atomic relaxations, 
while keeping the in-plane lattice parameters fixed to 
the relaxed bulk values. The supercell is large enough 
to accommodate the octahedral tiltings, and commensu- 
rate with the orthorhombic phase in the ab plane. The 
surface breaks the bulk rhombohedral symmetry, most 
importantly the three-fold axes of the i?3c phase, there- 
fore falling back onto the Prima tilt system. Doping is 
then introduced explicitly by substituting La atoms by 
Sr, with the resulting structure having a quasi-random 
distribution of 13 Sr atoms and 27 La atoms. The system 
is thick enough to ensure that the electronic structure at 
the center of the slab is similar to that of the bulk. 

We find that relaxation of the atomic positions gives 
rise to a buckling in the surface layer, with the Mn atoms 
shifted out of the plane. There is also a change in the oc- 
tahedral tiltings and of the Mn-0 bonding distances. In 
Fig. da) we plot the tilting angle of the MnC>6 octahedra 
as a function of the position of the Mn layer across the 
slab. It can be seen that in the central layers, the tilting 
approaches the value obtained for the relaxed bulk struc- 
ture (dashed line) while it reduces towards the surface. 
In addition, we find that the buckling of the Mn atoms 
with respect to oxygens in the surface layer propagates 
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FIG. 1: (a) Octahedral tilt angle versus position across the 
slab, z, in units of the cubic lattice constant, ao. The dashed 
line shows the value for the theoretically relaxed bulk geome- 
try, (b) Off-centering (in A) of the Mn atoms from the plane 
defined by the neighbouring O atoms in the Mn02 layer. The 
z origin is defined at the center of the slab. Error bars indicate 
the differences within the same atomic layer. 



to the subsurface giving rise to an off-centering similar 
to that of ferroelectric perovskites. The magnitude of 
this off-centering, shown in Fig. [TJb), decays monotoni- 
cally towards the center of the slab, where it is zero due 
to the cancelling effect of both surfaces. The distortion 
results in important changes in the Mn-0 bonding dis- 
tances along the z direction, whereas the in-plane bonds 
remain essentially unchanged. 

Although the Mn-0 buckling at the surface layer is 
a result of local surface chemistry (see the surface-state 
shift below), the off-centering deeper into the bulk is not. 
The situation is better described as a polar, ferroelastic 
instability, which is triggered by the surface buckling. 
Since this instability is not dominant in the bulk phase, 
it therefore decays from the surface into the bulk. Such 
surface-triggered decaying ferroelastic distortions are de- 
scribed in the literature [l9(, and can be very slowly de- 
caying 20]. Our calculations, however, allow us to as- 
certain only that the decay length A D is not smaller than 
~10A, since the antisymmetric character of the order pa- 
rameter [Fig. []Jb)] forces it to go through zero at the 
center of the slab. 

The length-scale of the decay towards the surface of the 
octahedral tilt order parameter, At, should be captured 
in our calculations, since the symmetry of the slab does 
not constrain it in this case. We find that it is of the or- 
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FIG. 2: (Color online): Partial density of states for up/down 
spins (top/bottom) for the Mn 3d orbitals in the center of 
the slab (a) , and the surface layer before (b) and after (c) the 
structural relaxation. The surface state is mainly composed 
of d z 2 and d xy orbitals, which are shown as shaded regions. 



der of a few atomic layers. However, we can not rule out 
the possibility that both order parameters are coupled at 
the surface, whereby A t would be in fact determined by 
the constraint on A Q . Indeed, couplings and/or competi- 
tions between polar ferroelectric modes and tiltings and 
rotations of oxygen octahedra are well known to occur 
in related perovskites 21, 22, 23| Interestingly, however, 
the decay of the untilting seems oscillatory and commen- 
surate with the crystal; both monotonic and oscillatory 
behaviors can be expected in surface ferroelastic decays 

Although not the main focus of this work, next we 
briefly compare the bulk and surface electronic struc- 
tures within the GGA. Fig. [2] shows our calculated Mn 
3d partial densities of states (PDOS) in the center of the 
slab (a), and for unrelaxed (b) and relaxed (c) surface 
configurations. At the center of the slab, the system is 
half- metallic, with e g character for the states at £p, and 
localized t\ g states around 1.5 eV below, in fair agree- 
ment with previous calculations [J] and photoemission 
experiments [24j ]. Importantly, this half-metallicity per- 



sists at the surface, and since the majority spin metallic- 
ity results from a broad Mn 3d - O 2p hybridized band 
(the Mn d z 2 and d x 2_ y 2 states, as well as the p z of the 
oxygen underneath the Mn, dominate the Fermi level at 
the surface) it is likely to be accurately represented in 
our GGA methodology. For the unrelaxed (cubic) struc- 
ture, the surface termination induces a localized surface 
state, of mainly d z 2 and d xy character. This state had 
been previously observed in LSDA and SIC-LSDA calcu- 
lations with unrelaxed cubic surface geometries, and orig- 
inates after the splitting of the states e g — » d z 2 1 d x 2_ y 2 
and t 
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d X y,(dyz,d X z) related to the surface-induced 
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symmetry breaking, from cubic Oh to C± v 
surface state shifts to lower energies away from ep when 
the geometry relaxes. This is the driving force for the 
surface buckling. 

Concerning magnetic properties, we observe that the 
zero kelvin magnetic moment remains almost unaltered 
up to the very surface, with a value of 3.55/xs at the 
relaxed surface compared with 3.53/is in the bulk. The 
ground state remains ferromagnetic. In order to address 
the finite temperature magnetism, we explored different 
magnetic configurations and found that the energetic cost 
of reversing one magnetic moment at the surface is re- 
duced to 30% of the bulk value, indicative of lower effec- 
tive magnetic interactions, and in line with experimental 
results [2(. Forces on the atoms of up to 1.54 eV/A (1.01 
eV/A ) have been found upon single spin reversal in the 
surface (in the bulk), indicating (i) that the surface ferro- 
magnetic coupling would weaken even more as compared 
to the bulk, and (ii) a possible important magnon-phonon 
coupling. 

Two competing considerations arise when relating the 
magnetic results to the structural ones. The double- 
exchange (DEX) mechanism, in which localized 3d ti g 
electron spins interact ferromagnetically through the 
e g conduction-electron hopping between nearest neigh- 
bours, is frequently used to describe the electronic prop- 
erties of manganites. On one hand, the electron-hopping 
energy among Mn e g states in this DEX model is ex- 
pected to decrease with the bending of the Mn-O-Mn 
bonds, in agreement with the observed reduction of the 
ferromagnetic Curie temperature Tm in doped mangan- 
ites [27(| . The octahedral untilting towards the surface 
discussed before would thus suggest a surface enhance- 
ment of hopping and Tm, at odds with the mentioned 
observations. On the other hand, however, the hopping 
energy should be much more affected by the substantial 
off-centering, thus altering the inter-layer coupling. 

The surface-triggered ferrodistortive instability dis- 
cussed above is very similar in kind and magnitude to 
that of ferroelectrics like BaTi03. It is thus tempting to 
consider multiferroic possibilities for this system. Con- 
sidering the calculated values for the atomic displace- 
ments, and taking ~ +3e for the nominal charge of Mn, 
we would obtain a polarization of ^10/iC/cm 2 at the sur- 
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face layer, comparable to the values obtained in similar 
Mn-based multiferroics. It must be remembered, how- 
ever, that it is a metallic system, and thus any polariza- 
tion will be screened by the itinerant carriers. In addition 
to the ferroelastic characteristics discussed above, the 
surface buckling could display surface piezoelectric sig- 
natures, considering the Thomas-Fermi screening length 
of this (poor) metal, affecting the work function and its 
strain dependence. This could be a relevant consideration 
for understanding the orientation of dipolar molecules de- 
posited on the surface. Spin valve devices based in LSMO 
have already been demonstrated with Alq3 (8-hydroxy- 
quinoline aluminum) (28l |. and a substantial shift in its 
molecular energy levels has been measured [29(. The 
LSMO surface dipole would favor the orientation of the 
Alq3 molecular dipole towards the vacuum side, and the 
resulting electrostatic potential would lower the molecu- 
lar levels as observed [29j . 

It should be remembered that the presented results 
here are for the defect-free Mn02-terminated LSMO sur- 
face. They may be very different for (La,Sr)0 termina- 
tion [26| . and, most importantly, in the presence of de- 
fects like oxygen vacancies. Recent photoemission exper- 
iments have suggested the presence of strongly localized 
Mn +2 states at the surface of LSMO, related to oxygen 
vacancies [lo| • DFT results on the effect of these defects 
in this LSMO surface will be discussed elsewhere [31 1 . 

In conclusion, our first-principles structural relax- 
ation of the LSMO (001) surface has revealed a sur- 
prising surface-triggered ferroelastic instability of an off- 
centering ferrodistortive kind. The distortion is driven 
by the surface electronic structure, and causes a weaken- 
ing of the tendency to magnetic order close to the sur- 
face. We hope that our observations will stimulate fur- 
ther work towards the characterization of the discussed 
instabilities, as well as possible exploitation of the newly 
observed surface polarization and its coupling to the mag- 
netic behavior. 
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